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ABSTRACT 



We present broad band spectra of a sample of 21 low luminosity sources in the 
Trapezium Cluster, with masses in the range 0.008 - 0.10 M© (assuming an age of 
1 Myr). These were selected for low extinction in most cases and are located west 
of the brighter nebulosity. The spectra are in the H bandpass (1.4-1.95 fim) and 
K bandpass (1.9-2.5 /itm) also for most of the brighter sources, with a resolution of 
50 nm. They were taken with the United Kingdom Infrared Telescope (UKIRT) us- 
ing the CGS4 spectrometer. Absorption by water vapour bands is detected in all the 
substellar candidates except one, which is a highly reddened object with strong H2 
emission and an anomalously blue (I- J) colour, implying that it is a very young cluster 
member with circumstellar matter. The observation of prominent water vapour bands 
confirms the low Effective Temperatures implied by our (I- J) colour measurements in 
an earlier paper and would imply late M or L spectral types if these were older field 
dwarfs. However, the profiles of the H bandpass spectra are very different from those 
of field dwarfs with similar water absorption strength, demonstrating that they are not 
foreground or background objects. In addition, the CO absorption bands at 2.3 /im 
and the Nal absorption feature at 2.21 fxm are very weak for such cool sources. All 
these features are quite well reproduced by the AMES-Dusty-1999 model atmospheres 
of Allard et al. (2000,2001), and arise from the much lower gravities predicted for the 
Trapezium sources (3.5 < log(g) < 4.0) compared to evolved objects (log g ~ 5.5), 
This represents a new proof of the substellar status of our sources, independent of the 
statistical arguments for low contamination, which are reexamined here. The very late 
spectral types of the planetary mass objects and very low mass brown dwarfs demon- 
strate that they are cluster members, since they are too luminous to be field dwarfs 
in the background. We also present additional UKIRT photometry of a small region 
in the south of the Trapezium cluster where the extinction and nebular brightness are 
low, which permitted the detection of objects with 1 Myr masses slightly lower than 
our previous least massive source at 8 M.j up . Following a minor update to our previous 
J band photometry, due to a new UKIRT filter calibration, there are ~ 15 planetary 
mass candidates in the full dataset. 

Key words: stars:low-mass, brown dwarfs, stars: pre-main sequence, stars: formation, 
(stars:) circumstellar matter 



1 INTRODUCTION 

Brown dwarf candidates have been discovered in large num- 
bers in many star formation regions and very young stellar 
clusters, through deep infrared photometric searches (eg. 
Comeron et al. 1993,1996; McCaughrean et al. 1995; Be- 
jar, Zapatero-Osorio & Rebolo 1999). Some of these surveys 
have yielded objects of planetary mass, in the Chamaeleon I 



molecular cloud (Tamura et al.1998; Oasa et al.1999) in the 
cr-Orionis cluster (Zapatero-Osorio et al.2000) and in the 
Trapezium Cluster in Orion (Lucas & Roche 2000, hereafter 
Paper I; Hillenbrand & Carpenter 2000). In most cases there 
is a very high probability that photometric candidates are in- 
deed substellar cluster members. Some candidates with near 
infrared JHK photometry exhibit K bandpass excesses due 
to hot circumstellar dust, which demonstrate their youth 
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(eg. Oasa et al.1999), although this criterion will miss the 
majority of young stellar objects (YSOs). In the Trapez- 
ium cluster, even substellar candidates without K band- 
pass excesses can generally be assumed to be cluster mem- 
bers because any measured extinction precludes confusion 
with foreground field dwarfs while the dense backdrop of 
the OMC-1 cloud and the high galactic latitude (b=18°) 
both mitigate against background contamination. 

It is still important to observe brown dwarf candidates 
spectroscopically, both to confirm their status and to ex- 
plore the rich physics of their atmospheres, which have much 
lower gravities than their evolved counterparts. There are 
effects specific to very young sources which complicate the 
interpretation of their spectra. These are: (1) light scattered 
and emitted by circumstellar matter, which modifies the 
signal from the photosphere; (2) reddening by the cluster 
medium, which can only be imperfectly removed; (3) the 
lower gravity of the atmosphere. The low gravities mean 
that the spectra cannot be simply compared to those of 
Class V dwarfs and such spectra have not been success- 
fully modelled in the near infrared until now. However, the 
new AMES-Dusty-1999 models of Allard et al.(2000, 2001) 
which employ the NASA Ames water line list and a more 
detailed treatment of dust, appear to be considerably more 
successful than previous models (see Section 4) . Evidence for 
low gravity spectral features has been seen in the Pleiades 
(Martin, Rebolo & Zapatero-Osorio 1996) which have an 
age of ~ 100 Myr. Luhman & Rieke (1998) and Luhman et 
al. (1998a) published optical spectroscopy of large samples of 
YSOs in IC348 and L1495E, including some substellar ob- 
jects. They allowed for the effects of low gravity by fitting 
spectra to templates formed from the average of Class V 
dwarfs and Class III giants, which produced a good match 
to the data. However this approach is not appropriate for 
the very low mass sources in the Trapezium and cr-Orionis, 
since the few good quality infrared spectra of late-M type 
Class III and Class I giants in the literature (eg. Lancon 
& Rocca-Volmerange 1992) do not have similar profiles to 
those presented here, and correspond to even lower gravity 
objects (—2 < log(g) < 1. Moreover there are no L-type 
giants to compare with the cooler sources found in Orion. 

Wilking, Greene & Meyer (1999) analysed K bandpass 
spectra of brown dwarfs in p Ophiuchus by comparison with 
Class V dwarfs, which may lead to errors in the derived effec- 
tive temperatures. We experienced the same problem in Pa- 
per I when deriving temperatures and luminosities by com- 
parison with Class V IJH colours. Fortuitously this does not 
appear to have caused significant errors because: (a) many 
sources with low extinction were observed with (J-H)— 0.5 
to 0.6, indicating that this colour does not increase to Red 
Giant Branch values ((J-H) ~ 0.9) in very young brown 
dwarfs and (b) the AMES-Dusty-1999 models indicate that 
low gravity causes only small changes in the broad band 
colours, comparable to our typical photometric uncertainty 
of 0.1 magnitudes. 

In this paper we present low resolution H and K band- 
pass spectra of a sample of Trapezium cluster members in- 
cluding low mass stars, brown dwarfs and planetary mass ob- 
jects. The spectra are analysed by comparison with both the 
AMES-Dusty-1999 models and published spectra of Class V 
dwarfs. Optical spectra will be presented in a future paper. 
New photometric data are also presented covering the south- 



ern region of the Trapezium cluster. We also readdress the 
arguments concerning the nature of these sources, in light 
of the spectroscopic data. 

1.1 Note on Nomencalature 

In this paper we use the less controversial but not very de- 
scriptive term 'planetary mass objects' to describe sources 
below the deuterium burning threshold. Some in the extra- 
solar planets community object to the term 'free floating 
planets' to describe objects which are not orbiting a star 
and probably form in a molecular cloud core in a manner 
similar to stars (although the word planet originally meant 
'wanderer'). As a possible alternative, we suggest the more 
compact term 'planetar' as a contraction of 'planetary mass 
object believed to have formed like a star'. A large amount 
of theoretical work is now in progress, examining a variety of 
different formation scenarios. There may prove to be more 
than one mode of formation for these low mass bodies, as is 
believed to be the case for stars. 



2 OBSERVATIONS 
2.1 Spectroscopy 

Long slit spectroscopy of sources in the Trapezium clus- 
ter was carried out at the United Kingdom Infrared Tele- 
scope (UKIRT) on 27-30 November 1999, using CGS4, a 
near infrared Cooled Grating Spectrograph. The observer 
was PWL. Sources were observed in pairs aligned along the 
slit of the spectrograph and the sample was selected from 
the dataset of Paper I and the new imaging data presented 
here, with the intention of sampling the mass range from 
0.008 to 0.10 Mq, i.e. from planetary mass to very low mass 
stars. 21 such sources were observed and 2 sources of higher 
mass. There were the following additional criteria. (1) Lo- 
cation in a region of low nebular surface brightness in order 
to improve sensitivity. (2) Low extinction towards at least 
one member of a pair, defined by (J-H) < 1.5 or A(V) ^ 7.5, 
also to improve sensitivity. (3) Most sources were selected for 
'normal' (I- J) colours in an attempt to minimise contamina- 
tion by light scattered from circumstellar matter (see Paper 
I) but 2 sources with the anomalously blue colours shared by 
14% of the cluster population were observed to see whether 
their spectra were different from the others. (4) Location 
within ~ 40" of a sufficiently bright star (J mag < 15) to 
permit precise location of optically invisible targets within 
the slit by offsetting to the target. In practice, all 21 of the 
low mass targets were located at least 105 arcsec west of 6i 
Orionis C at the cluster centre. 

CGS4 has a 256x256 InSb array sensitive from 1-5.5 /im 
and a spatial scale of 0.62 arcsec/pixel. The 40 1/mm grating 
was used in first order, which yields a spectral dispersion of 
25 nm/pixel but a 1.2"wide slit was employed in most cases, 
yielding an instrumental profile with a width of 50 nm. This 
corresponds to a spectral resolution of R=330 in the H band- 
pass and R=440 at K. The wide slit was used to ensure that 
most of the flux from both members of a widely separated 
pair would fall in the slit. Each spectrum covers a 0.62 ^m 
wavelength range. The H bandpass spectra were centred at 
1.65 /im and the K bandpass spectra at 2.2 /xm, overlapping 



© 2001 RAS, MNRAS 000, |-?? 



Infrared Spectroscopy of Substellar Objects in Orion 3 



[atm228nl]: BS 1687 INDEF ap:l beam:l 




Wavelength (microns) 

Figure 1. Relative transmission of Atmosphere and Spectrograph 
per unit wavelength. 

in the 1.9 fim region. One pair of sources was observed at 
25 nm resolution with a 0.6"slit, though not in very good 
conditions; no additional narrow line features were revealed. 
This pair of sources was observed with H bandpass spectra 
centred at 1.58 /im on the first night, in attempt to include 
the steep J band water edge at 1.34 /j,m, but the instrument 
efficiency in this configuration was found to be essentially 
zero for wavelengths A < 1.4 /im, so the longer wavelength 
setting was used thereafter. The relative throughput of the 
instrument and the atmosphere in the H bandpass is shown 
in Figure 1. The observations and source notes are detailed 
in Table 1. For all other pairs, the array was stepped over 2 
pixels in a whole pixel step to remove bad pixels and the 
telescope was nodded along the slit in order to subtract 
the background, with background-limited individual inte- 
grations of 30 to 60 s duration. Total exposures time ranged 
from 10 minutes for the brightest sources to 1.5 hours for the 
faintest. 21 primary sources were observed in the H band- 
pass and 11 of these in the K bandpass also. We also present 
H bandpass spectra of 2 fairly low mass sources which ap- 
peared in the long slit by chance. The smaller number of 
K bandpass observations was due to constraints on observ- 
ing time (weather and computer failures) . Priority was given 
to H bandpass observations because of anticipated compli- 
cations owing to emission by hot dust in the K bandpass 
and because the water absorption troughs are slightly more 
prominent in the H bandpass. A small number of red main 
sequence stars were also observed in order to check the reli- 
ability of the data. Their spectra agree very well with those 
published by Leggett et al.(2000) for the same objects, also 
obtained with CGS4 at UK1RT. 

The data were coadded into background subtracted 're- 
duced groups' using the CGS4dr software package and spec- 
tra were optimally extracted using the TWODSPEC package 
in IRAF, the Image Reduction and Analysis Facility. Argon 
arc spectra, tungsten lamp flatfields and standard stars were 
observed frequently to compensate for changes in the wave- 
length calibration when the slit orientation was changed and 
to accurately correct for the changing telluric water vapour 
column density. F-type stellar standards were generally used 
in preference to A-types since they have weaker HI absorp- 
tion lines. The fairly weak Brackett series lines were removed 
from the stellar standards by fitting the continuum with 
a high order cubic spline, dividing the fit by the data to 



produce a line spectrum, and multiplying the data by the 
line spectrum. The strong Paschen a line at 1.875 /xm co- 
incides with a deep trough due to telluric water absorption 
and could not be adequately removed, so we have removed 
a small section of the spectra near this wavelength. The 
Brackett series limit occurs at 1.46 fim and coincides with a 
narrow telluric feature at the short wavelength edge of the 
H bandpass, making it impossible to accurately remove. In- 
spection of the NextGen models of Hauschildt et al.(1997) 
indicates that the Brackett series limit is weaker than the 
telluric feature in F-type standards and it does not appear 
to be noticeable in the spectra presented here. 

In addition to background shot-noise, a significant 
source of noise arises from the spatial gradients in the nebu- 
lar surface brightness, which make the background subtrac- 
tion difficult. Owing to these gradients the reduced group 
images are not flat, but show the effects of subtracting two 
frames separated by distance of the nod along the sky: either 
10 or 20 arcsec was used depending on the separation be- 
tween the members of a pair. The images also show the spec- 
tral residuals due to fluctuations in the strength of telluric 
OH emission lines. A polynomial was interactively fitted to 
the residual background at each wavelength and subtracted, 
using the APALL package in TWODSPEC, and this proce- 
dure was generally very successful in removing both nebu- 
losity gradients and residual OH emission. The problem was 
greatest for the fainter sources, whose peak brightness is 
comparable to fluctuations in nebular suface brightness on 
a scale of a few arcsec, even for sources selected in regions of 
low nebular background. Where necessary, 2 iterations were 
used to produce a flat background. For the faintest sources, 
the background may be slightly over or under subtracted, 
but agreement between spectra taken at the two nod po- 
sitions was used as a check to prevent this occuring. Only 
one spectrum was lost due to this problem: the K bandpass 
spectrum of the faintest source, 051-147, was rendered use- 
less by a small step in the surface brightness on the same 
scale as the nod length, but an H bandpass spectrum was 
successfully extracted. We recommend that a nod length of 5 
arcsec or less be used in star formation regions when taking 
long slit spectra of faint sources. 

2.2 Photometry 

Additional photometry in the IJH bands was obtained at 
UKIRT, using the infrared camera UFTI, in the manner de- 
scribed in Paper I, using 900 s exposures. The new data were 
taken on the nights of 11-12 November 1999 by observatory 
staff. They consisted of three contiguous 1.5 arcmin UFTI 
fields located 3 arcmin south of 6\ Ori C, oriented east-west 
as shown in Figure 2. This region was selected because our 
existing data indicated that the nebular surface brightness 
and average extinction would be low in this region, and pub- 
lished K bandpass maps showed that the stellar density is 
reasonably high (Hillenbrand 1997). These features assisted 
in the detection of some new very faint sources, as described 
in Section 3.1. The seeing conditions were fairly good, with 
a FWHM of 0.5-0. 6"in all three filters, aided by the tip/tilt 
image stabilisation of the UKIRT secondary mirror. 

Photometry was performed using the DAOPHOT 
crowded field photometry package in IRAF as in Paper I, us- 
ing both manual photometry and automated crowded field 
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Figure 2: J Bandpass image of the full photometric dataset. The region of new data is below the horizontal line. The spectroscopic 
sample was taken from sources to the west of the brightest nebulosity. The locations of very low mass candidates (dereddened mag 

J,j r > 17.5, or M^S 0.016 Mq) are marked with circles. 



photometry. The former was generally preferred in order 
to accurately subtract the highly structured nebulous back- 
ground, except when the image profiles of sources overlapped 
significantly. Several sources were contained in the overlap 
region between the new data and the dataset of Paper I. The 
measured magnitudes typically differed by 0.1 mag or less, 
which is approximately the photometric error described in 
Paper I. The main source of error in the brighter sources is 
the image profiles, which are spatially variable in the over- 
lap region since this is t the very edge of both the old and 
new datasets. Each UFTI field is made up of a 9-point image 
mosaic with small jitters. 

2.2.1 Calibration Update 

The J band magnitude system for UFTI at UKIRT has re- 
cently been updated by observatory staff, who discovered 
that the transformations from the old IRCAM system to the 
new UFTI system are different for cool red stars than other 
main sequence stars. This is because these stars have water 
absorption bands in similar locations to the water bands in 
the earth's atmosphere and UFTI employs the new J band 
filter developed for the Mauna Kea Consortium, optimised 
for the low water column above this observing site. The new 
transformation has been thoroughly calibrated via measure- 
ments of late M type and L type stars. Since our Paper I 
dataset was calibrated using red stars (in order to calculate 



transformations from Cousins I to the UFTI I band) this 
requires us to update our J band photometry, and conse- 
quently the I band transformations also (see UKIRT web 
pages for the details). The new J band magnitudes are typ- 
ically 0.1 magnitudes brighter than before, which leads to 
a slightly lower luminosity for sources with significant ex- 
tinction, since the calculated extinction is reduced. Sources 
with near zero extinction are slightly more luminous than 
had been thought. However changes of this order are of very 
little consequence, since the mass-luminosity relation is quite 
steep throughout the brown dwarf regime and the change is 
similar to both the typical photometric uncertainty and the 
theoretical uncertainty in the mass-luminosity conversion. 

The bolometric correction, BCj, formerly used to con- 
vert dereddened J magnitudes, Jdr, to luminosities is also 
updated in consequence: 

BCj = 2.01; J dr > 13.71 (la) 
.BCj = 0.209Jdr - 0.855; J dr < 13.71 (16) 

The new models used here predict the fluxes in each 
filter directly without the need for a bolometric correction. 



3 RESULTS 
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3.1 Photometry 

43 new sources were detected at both J and H bands in the 
new data, of which 31 were detected at / band. A colour- 
luminosity plot, with the new data added to that of Paper 
I, is shown in Figure 3 for all 557 unsaturated sources. We 
plot 2 different 1 Myr isochrones as almost vertical solid 
lines at the left of the picture, to which sources may be 
dereddened. Both of these are self consistent predictions us- 
ing the same model atmospheres for colour predictions and 
to provide the boundary condition for the evolutionary cal- 
culations. This removes the need to use uncertain bolomet- 
ric corrections. Furthest to the left is an isochrone supplied 
by Baraffe et al.(2001, in preparation), using the non-gray 
AMES-Dusty-1999 model atmospheres (Allard et al.2000, 
2001). The isochrone slightly to the right is the prediction 
of Baraffe et al.(1998), recently extended to 1 Myr and to 
the deuterium burning threshold (Chabrier et al.2000). This 
second isochrone has evolution and colours based upon the 
NextGen model atmospheres. It is useful to see both be- 
cause the Nextgen atmospheres correctly predict the colours 
of main sequence stars with effective temperatures above 
2500 K, whereas the AMES-Dusty models are successfully fit 
the colours of field brown dwarfs, below 2200 K, where dust 
has an important influence on the photosphere. The tem- 
perature at the deuterium burning threshold lies near the 
low temperature end of the NextGen regime but we cannot 
be sure which model will predict the (J-H) colours better 
in very young, low gravity photospheres. The true colour 
probably lies between the 2 predictions, as discussed below. 

Recalibrated photometry of the full dataset is available 
electronically^. The O'Dell & Wong (1996) coordinate based 
naming system is used as in Paper I. The advantageous fea- 
tures of the new dataset, as described above, permitted the 
detection of some very low luminosity sources: 4 new plane- 
tary mass candidates, 3 of which appear to be less massive 
than any of the sources described in Paper I, assuming an 
age of 1 Myr and employing either the isochrones of Bur- 
rows et al.(1997, hereafter BM97) or the new isochrone of 
Baraffe et al. However only 1 of these 3 faint sources has 
good photometry: 189-659 has m.j = 19.23, (J ' - H) = 0.63, 
which indicates a notional mass of 6.2 Mj up at 1 Myr (on 
the BM97 isochrone) and apparently negligible extinction. 
The other 2 sources have large photometric uncertainties in 
both the J and H bands, so their masses are highly uncer- 
tain. Intriguingly, these 3 faint sources are all located close 
together on the sky, in a 1 arcmin square region 3 arcmin 
due south of 9i Ori C, on the far side of the Orion Bar from 
the cluster centre. The Orion Bar is the ionisation front of 
the OMC-1 cloud, viewed edge-on as it curves into the line 
of sight. Sources which are located beyond it lie on lines of 
sight passing deep into the molecular cloud. Hence we might 
expect such sources to be younger than the rest of the clus- 
ter and hence less massive for a given luminosity. However, 
the very low extinction toward 189-659 indicates that this 
source probably lies in front of the cloud and might possibly 
be a minimum mass star in the foreground. The locations 
of all the sources with dereddened J magnitudes Jd r > 17.5 



* Photometry and spectra are available via anonymous FTP to 
star.herts.ac.uk, in pub/Lucas/Orion. 







1 1 1 1 

1 


1 






'.■ ; 
7 • • 


A(V) = 7 




• 


T^S* •.« • • 

.»*• h .• • . ... • . . 

•»\,\ t • < . 
•r ••• **\ • • • * • • " 

• >\- . . . .-. •• • 

K*bv'- -\- •• • 


•• 






? ss.* ° v\ ° 

V* ° \\f 

\> \\ 
\f \\ 


1 



2 

(J-H) 



Figure 3: Colour-magnitude diagram for the 557 unsaturated J 
bandpass sources. Two theoretical 1 Myr isochrones are plotted at 
left, with reddening lines indicated for 0.08 Mq and the 0.012 Mq 
deuterium threshold. Open circles represent sources with large J 
band photometric uncertainties (>0.2 mag). 



(i.e. M< 0.012 to 0.016 M Q at 1 Myr) are plotted in Figure 
2. 

There are now approximately 15 planetary mass can- 
didates in the dataset, following the filter recalibration. The 
number of sources below the deuterium burning thresh- 
old cannot be specified very precisely, mainly because 
of the uncertainty in the (J-H) colour of their photo- 
spheres mentioned earlier. Most evolutionary calculations 
are in very good agreement concerning the minimum mass 
for deuterium burning, at 0.012-0.013 Mq (Saumon et 
al.1996; BM97; Chabrier et al.2000). In addition, the 1 Myr 
isochrones of BM97, D'Antona & Mazzitelli (1998) (herafter 
DM98, update to 1997 publication) and Baraffe et al.(1998) 
are in fair agreement (within 0.2 to 0.3 mag) with regard to 
source luminosities for low mass brown dwarfs. However, the 
small uncertainty in the (J-H) colours can slightly change 
the number of planetary mass objects by moving the lines 
of constant mass in Figure 3 to the right or to the left 
by 0.1-0.15 mag, since a number of sources appear to lie 
very close to the deuterium burning threshold. Applying the 
new Baraffe et al. isochrone, there are ~13 planetary mass 
objects (below the M< 0.012 Mq threshold calculated by 
Chabrier et al.), whereas the NextGen-based Baraffe et al. 
ischrone indicates there are 16 such objects. Using the main 
sequence colours and bolometric corrections as in Paper I, 
and the BM97 isochrone for luminosity, indicates that there 
are 15 such objects. Given that the J bandpass photometric 
errors are typically 0.1-0.2 mag for these faint candidates, 
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the important issue is not the precise number but whether 
the IMF extends below the deuterium burning threshold, 
which appears to be the case. 

The NextGen colour predictions, using the Miller- 
Tennyson water line list, are only very slightly too red for 
main sequence stars but these atmospheres poorly repro- 
duce the observed profiles of both field dwarfs with strong 
water absorption (Leggett et al. 2001) and the low gravity 
spectra presented here. The AMES-Dusty-1999 colour pre- 
dictions, using the NASA AMES water line list (Partridge 
& Schwenke 1997), are 0.1 to 0.2 mag too blue in main se- 
quence stars (see Allard, Hauschildt & Schwenke 2000) but 
reproduce the profiles of the spectra presented here quite 
well. A small colour discrepancy in (J-H) would appear to 
arise from incompleteness in the water line list in the 1.7- 
2.2 fj,m water absorption band, which is shown directly for 
Trapezium sources in Section 4.1. 

Note that we have made a conservative choice of red- 
dening law (see Paper I). The choice has little effect in any 
case because the planetary mass candidates necessarily have 
low extinction and because extinction laws have little varia- 
tion in the infrared (see Cardelli, Clayton & Mathis (1989)) 
except in extremely dense environments such as the accre- 
tion disks of embedded protostars, which may contain dust 
grains larger than 1 /j,m in radius. 

There may be small but unquantified effects on the (J- 
H) colours due to scattering by circumstellar matter and a 
small percentage of emission by hot circumstellar dust in 
some cases. Hot dust would make sources appear slightly 
redder (see Section 4) and overmassive; the direction of any 
scattering effect depends very much on assumptions about 
the circumstellar structure and viewing angle. The addi- 
tional uncertainties of cluster membership and candidate 
ages are discussed in Section 6. 

Erratum 

We note that 1 planetary mass candidate from Paper I 
was found to be double counted owing to a typing error. An 
internal cross check of all the coordinates shows that this is 
the only such case in the dataset of 557 unsaturated sources 
detected in the J and H bandpasses. 



an exposure. The noise is most obvious in very faint sources 
such as 084-104 (Figure 4(c)) in which many spurious narrow 
features remain despite some smoothing. The data quality is 
poor in the 1.8-1.95 /im region owing to strong telluric wa- 
ter absorption and also becomes less reliable at A < 1.5 fim, 
due to both telluric water absorption and the declining ins- 
tument efficiency. Some of the spectra are strongly reddened, 
with the result that the water absorption appears weaker on 
the long wavelength side of the peak and stronger on the 
short wavelength side, eg. object 013-306. In most cases the 
effect is minor owing to our preferential selection of sources 
with low extinction. The extinction toward each source is 
calculated using the IJH photometry, by comparison with 
measured main sequence star colours as described in Pa- 
per I, and noted in Figure 4(a-c). The DEREDDEN task in 
IRAF was used, which is based on the Cardelli, Clayton & 
Mathis extinction law, with reddening parameter R^=5.5 
as before and A(A) oc A -1 ' 61 in the infrared. This is a some- 
what imprecise procedure, because our typical photometric 
uncertainty of 0.1 mag in each filter leads to an uncertainty 
of 0.14 mag in the (J-H) colour, which corresponds to an 
uncertainty of 0.25 mag in H band extinction and 1.2 mag 
in V band extinction. However, the errors will be less for the 
brighter sources in the region of the spectroscopic sample, 
which was selected for faint nebulous background. The use of 
main sequence colours gives similar results to the NextGen 
atmospheres: the atmospheric models predict changes of no 
more than ~ 0.1 mag in (J-H) due to lower gravity in young 
sources, so we do not think this uncertainty will be larger 
than measurement error. 



3.2 Spectroscopy 

3.2.1 H band 

21 H band spectra were taken of sources with masses M 
<O.1M0. These are plotted in Figure 4(a-c), as F„ spectra 
in order to suppress the apparent noise at short wavelengths. 
Dereddened spectra are plotted above the observed spectra 
for those sources with non-zero extinction. These spectra 
are not flux calibrated since for some observations observ- 
ing conditions were not photometric. The fainter sources 
are Gaussian-smoothed to enhance low resolution features. 
Nearly all of the spectra are dominated by strong water ab- 
sorption bands to either side of a sharp peak located between 
1.68 and 1.70 /im. The exceptions are plotted in the lower 
part of Figure 4(b). No real narrow spectral features are 
unambiguously detected in the H bandpass. Pixel to pixel 
variations therefore illustrate the noise in each spectrum, 
which has wavelength dependent contributions from narrow 
telluric OH emission lines and broad band telluric water ab- 
sorption, both of which vary on a timescale of minutes during 
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Figure 4(a): Observed bandpass spectra of Trapezium sources, with dcreddencd spectra overplotted as the upper line for each source 
when extinction is non-zero. Masses derived from the BM97 1 Myr isochrone are indicated. The noise is indicated by the pixel to pixel 
variations, since no narrow spectral features are detected in the H bandpass. 
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Figure 4(b): Observed and Dereddened H bandpass spectra, as Figure 4(a). 047-436 shows no water absorption and 2 others have very 
weak water absorption. 
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Figure 4(c): Observed and Dereddened H bandpass spectra for the least massive candidates, with some gaussian smoothing applied. The 
apparent absence of flux at A< 1.5 ^m in 023-115 is attributed to the decline in signal/noise at those wavelengths, due to atmosphere, 
instrument efficiency, extinction and (presumed) water absorption in the source itself. 
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The peaked H bandpass profiles in Figure 4(a-c) are 
very different from those of late M and L-type dwarfs in the 
Local Neighborhood with strong water absorption bands, as 
shown by our spectra of 2 red standards in Figure 5(a-b), 
compared with a composite spectrum of the Orion sources. 
The composite spectrum is an average of 12 reasonably high 
quality spectra with fairly strong water absorption, weighted 
according to the apparent noise in each spectrum. The com- 
posite has slightly weaker water absorption than the average 
of the sample, since it is weighted toward warmer, brighter 
sources near the stellar/substellar boundary. Comparison 
with the data of Leggett et al.(2001) (overplotted in Fig- 
ure 5(a)) for the same red standards demonstrates that the 
spectra are reliable and the differences from field stars are 
real. In spectra of field dwarfs there is strong water absorp- 
tion at A < 1.55 fim and A > 1.70 /im, but a shallower slope 
in between, which appears as a flat plateau in Fa plots of 
the data (Figure 5(b)). This is clearly seen in the Fa plots of 
a large sample of late M and L dwarfs presented by Leggett 
et al.2001. marked by clear changes in slope at 1.70 /im 
and in the 1.55-1.60 fj,m interval. The peaked profiles of the 
Trapezium objects, which show no plateau, appear to be a 
consistent signature of youth and cluster membership. They 
are modelled in Section 4. There are no other clearly seen 
atomic or molecular features in the Trapezium sources at H 
band, which is unsurprising given the low spectral resolu- 
tion and modest signal to noise. Atomic absorption lines of 
Na and Ca have been observed in field dwarfs, eg. Jones et 
al.(1994) in high signal to noise spectra but many atomic 
lines appear to be weaker in very young brown dwarfs (eg. 
Bejar, Zapatero-Osorio & Rebolo 1999.) The water absorp- 
tion depths are quantified in Table 2, and discussed in Sec- 
tion 4. 

Only 3/21 sources have weak or zero water absorption. 
019-108 and 066-433 suffer very high extinction (A(V) = 
8.6, 12.6 respectively) and have very weak water absorp- 
tion, consistent with late K-type or very early M-type spec- 
tra. These earlier spectral types are not consistent with the 
temperatures expected at age 1 Myr even for these rela- 
tively luminous sources. (At 1 Myr masses, M=0.090 Mq, 
0.066 M and T e// =3400 K, 2800 K respectively are de- 
rived photometrically, using the BM97 isochrone). There- 
fore the H band spectra indicate that they are either even 
younger and hotter systems (age ~ 0.1 Myr) or they are 
background field stars. The dereddened fluxes of these stars, 
are faint enough to be consistent with late K-type dwarfs 
at d> 440 pc (mj = 13.72, 14.45 for 019-108 and 066-408 
respectively) judging by comparison with both the Lumi- 
nosity Class V sequence of Tokunaga (2000) and the models 
of BCAH for low mass stars at late K dwarf temperatures. 
019-108 would have to lie just behind the OMC-1 cloud at 
440 pc, while 066-433 would lie at only a slightly greater dis- 
tance. It therefore seems slightly more probable that these 
2 hotter, highly reddened sources are background stars but 
the alternative of extreme youth is not excluded, especially 
since the stellar mass proplyd source 178-232 shows a similar 
spectrum (see below). 

Only one source, 047-436, shows no water absorption 
in the H bandpass, simply a very red continuum in Figure 
4(b). This source is one of the 14% from Paper I with an 
anomalously blue (I- J) colour (though not an extreme case, 
1=18.277, (I-J)=3.21) and so its extinction cannot be well 



measured. At present, we attribute such colours to light scat- 
tered by circumstellar matter. An attempt to apply the stan- 
dard procedure leads to an insufficiently dereddened spec- 
trum, shown in Figure 4(b). The youth and cluster member- 
ship of 047-436 is demonstrated by the K bandpass spec- 
trum, discussed below. 

2 additional H bandpass spectra of more massive 
sources are presented in Figure 6. 1 of these, 017-410 is a 
YSO with a mass of ~ O.3M0 which appeared in the slit 
by chance, owing to the crowded nature of the cluster. It 
shows fairly strong water absorption but more on the long 
wavelength side of the 1.7 /J,m than the short wavelength 
side, which is different from the less massive sources, which 
may reflect the higher gravity of this source. Luhman et 
al.(2000) have obtained K bandpass and optical spectra of 
stellar mass Trapezium sources but there do not appear to 
be many H bandpass spectra of Trapezium stars in the lit- 
erature to use for comparison. The second source is 178-232, 
a more massive star (M~ O.4M0) with an anomalously blue 
(I- J) colour around which an emission line 'proplyd' struc- 
ture was resolved by the Hubble Space Telescope (O'Dell & 
Wong 1996). This has a very red spectrum, perhaps due to 
extinction by circumstellar matter rather than dust in the 
cluster medium and shows relatively weak water absorption, 
which would be expected for a hotter, younger object. 

3.2.2 K band 

11 K band spectra are shown in Figure 7(a-b), and dered- 
dened H+K spectra are plotted in Figure 8(a-b). The H 
and K spectra are joined by averaging the overlap region 
from 1.89-1.97 fiva, with some uncertainty (~ 10%) in the 
flux level owing to the telluric water absorption at those 
wavelengths. 9/11 of the K band spectra show strong water 
absorption bands on either side of a maximum near 2.25 /im, 
as expected for young brown dwarfs. The other 2 spectra, of 
019-108 and 047-436, are consistent with very weak or zero 
water absorption, similar to the H band data. The noise is 
indicated by the pixel to pixel variations as in the H band- 
pass data. 

The v=2-0 vibration-rotation bands of CO at A > 
2.29 fim are detected in 4/11 sources (see Table 3) via the 
edge at 2.29 /J,m and marginally detected in one other source. 
Only one source, 053-503, shows the multiple bandheads 
clearly at 2.29, 2.32, 2.35 and 2.38 /jm but they are also 
visible in 014-413. The weakness or absence of CO absorp- 
tion is surprising for the low temperature photospheres in- 
dicated by the strong water absorption. However, it appears 
to be a feature of the youth and low gravity of cluster mem- 
bers, as demonstrated by the AMES-Dusty-1999 models in 
Section 4. Veiling by hot circumstellar dust cannot explain 
it, since the spectra decline steeply toward 2.5 /im, show- 
ing that water absorption is not significantly veiled at these 
wavelengths. 

At K band we might expect that hot dust would redden 
the spectra but if this is occurring the effect is too small to 
be unambiguously identified, except in the cases of 014-413 
and 047-436, the only sources in Figure 8(a-b) showing much 
more flux in the K bandpass than at H . 014-413 appears to 
be associated with a weak wind collison shock front in the 
data of Bally, O'Dell & McCaughrean (2000) but shows no 
evidence of emission lines in the spectrum shown in Figure 
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Figure 5(a-b): Red Standards and Orion composite spectrum. (a)(left) F„ spectra. Leggett et al.(2001) data are overplotted as a solid 
line for the field dwarfs to demonstrate the reliability of the spectra, (b) (right) spectra plotted as F^, illustrating the plateau seen in 
field dwarfs in the H bandpass. The composite is a weighted average of the warmer, more massive sources, so it has slightly weaker 
water absorption than the average of the sample. Weak narrow absorption features in the composite are generally attributed to telluric 

OH lines. 



8(a). The possible effect of hot dust is shown in Figure 11, 
see Section 4. We do not have K band photometric data 
but Hillenbrand & Carpenter (hereafter HC) published HK 
photometry of the central region of the Trapezium Cluster. 
In combination with extinction values derived from our I JH 
data, these data should identify any sources with large K 
band excesses. The HC data contains photometry for only 2 
of the 11 sources in our K band spectroscopic sample (which 
lies to the west of the cluster core). Their H band fluxes of 
HC agree with our data to within 0.18 magnitudes for the 
2 sources, which is broadly consistent with the photometric 
uncertainties but may indicate T Tauri-like variability. Nei- 



ther of these has any K band photometric excess, which is 
consistent with the spectra. Some sources in the much larger 
photometric dataset do appear to have K bandpass excesses 
(see Section 6 and Table 4). Wilking, Greene & Meyer calcu- 
lated the contribution of hot dust using JHK photometry 
for low mass stars and brown dwarfs in the p Ophiuchus 
star forming region. They found that hot dust contributed 
0% to 20% of the total flux in the K bandpass in the major- 
ity of cases, which would have only a small effect on spectral 
profile. Star formation is nearly complete in the Trapezium 
Cluster and the average age of this sample of low-extinction 
sources is probably greater than that of sources in p Ophi- 
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Figure 6. Observed H bandpass spectra of 2 more massive sources (0.3-0.4 Mq). 178-232 is a proplyd; the weak water absorption and 

red spectrum is consistent with an embedded, warmer source. 



uchus so it is not too surprising that we do not see many 
large infrared excesses here. 

The Nal line at 2.21 /jm is the strongest atomic aborp- 
tion line in the H and K band spectra of field dwarfs (see 
Jones et al. 1994). It is clearly visible in the K band spectra 
(see Figure 16) but is not seen in the Trapezium sources. We 
would expect it to be seen in the high quality spectrum of 
053-503 (see Figure 16) but there is only a small local min- 
imum near that wavelength, which cannot be distinguished 
from noise. This weakness of this feature is also predicted by 
the AMES-Dusty-1999 models as a signature of low gravity, 
as described in Section 4.1.2. 

The K band spectrum of 047-436 (Figure 7(b)) has a 
very red continuum in which water absorption is weak or 
absent, following the trend seen at H band. However, there 
is a rich sequence of H2 emission lines, with strong lines at 
2.12, 2.41 and 1.96 /im and several weaker lines apparent 
also (see Table 2). As mentioned above, this is an anoma- 
lously blue source in the I band and these are usually also re- 
solved as proplyds by the Hubble Space Telescope (O'Dell & 
Wong 1996), as discussed in Paper I. This particular source 
is listed as being stellar in the data of O'Dell and Wong 
however, perhaps because it is located a large distance (3 
arcmin, or 0.4 pc) from the cluster centre, leading to less 
excitation of the optical emission lines that mark a proplyd, 
and insufficient contrast with the faint nebulous background 
for a silhouette disk to be observed. 047-436 lies < 1 arcmin 
north of the complex of bow shocks associated with LL Ori., 
observed in Ha by Bally et al.(2000). 

The H2 emission is the signature of a very young sys- 
tem with plentiful circumstellar gas. Such a signature can 
be caused by shocked gas in a molecular outflow or by UV 
excitation of gas in a circumstellar accretion disk or enve- 
lope i.e. a proplyd similar to those discovered by O'Dell & 
Wen (1994). In Orion the UV radiation is supplied mostly 
by 6>i Ori C. Chen et al.(1998) resolved the 2.12 /im H 2 from 
2 proplyds with NICMOS and found that it was confined to 
the subarcsecond disk silhouettes seen in continuum imag- 
ing. However, the emission from the 3 strongest H2 lines in 
047-436 is extended along the slit of the spectrograph by at 
least 3 arcsec (1300 au) in a southerly direction, the peak 
of the emission being coincident with the continuum source. 



There is no extension in the continuum and seeing condi- 
tions were fairly good (~ 0.6 arcsec) at the time the spec- 
trum was taken. The slit was oriented 18° west of north. The 
large extension of the H2 emission in 047-436, relative to the 
Chen et al. proplyds, may be due to the high extinction and 
much larger projected distance from the Trapezium stars, 
which would allow molecular material to survive to survive 
for a greater period of time. The reduced group image is not 
very sensitive to larger scale emission because a nod length 
of only 10 arcsec was used in taking the spectra and sub- 
tracting the background. However, the extended emission 
appears to be weaker in the northerly direction, extending 
for < 1.5 arcsec. The strength of the emission in the south- 
ern direction declines monotonically with distance from the 
peak, dropping by a factor of ~2 at a distance of 2 arcsec 
from the central source. There do not appear to be any H2 
emission line spectra of proplyds in the literature, although 
there are spectra of shocked gas in the vicinity of the BN 
object. Theoretical models of optical proplyd spectra, eg. 
Storzer & Hollenback (1999) predict that some of the lines, 
eg. the (1-0) S(l) line at 2.12 jj,m, may be collisionally ex- 
cited in the dense environment of the circumstellar disk, not 
excited by UV fluorescence. This makes it more difficult to 
understand the relative strengths of the emission lines and 
a full analysis is beyond the scope of this paper. However, 
the ratio of the (2-1) S(l) and (1-0) S(l) line strengths in 
Table 2 does suggest UV fluorescence. 

We conclude that 047-436, the only source entirely lack- 
ing H bandpass water absorption, is definitely a member of 
the Trapezium cluster. The slope of the 'dereddened' spec- 
trum in Figure 8(b) is not very reliable because it is an 
anomalously blue source (I- J) source. It is quite possible 
that hot dust contributes to the K band spectrum in this 
very young system but the bulk of the very strong redden- 
ing in the H bandpass must be due to extinction, perhaps 
by circumstellar dust. A very young system (age of order 
0.1 Myr) would have a hotter photosphere than the rest of 
the sample, so any water absorption would be weak and in 
this case veiled by extinction. The mass of 047-436 can be 
estimated as 0.06 Mq from IJH photometry in Paper I, as- 
suming an extinction of A(V)=5.5. A 0.1 Myr source would 
be less massive than a 1 Myr source for a given luminosity 
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Figure 7(a): Observed K bandpass spectra of Trapezium sources. The noise is indicated by pixel to pixel variations, as in Figure 4(a-c). 
069-210 and 055-231 are gaussian smoothed, because they were observed in poor conditions. 



© 2001 RAS, MNRAS 000, 



14 P.W. Lucas, P.F. Roche, France Allard, and Peter. H. Hauschildt 



. r tr ■ ii- ii mcr . F .i 



017-436 




S 



I I : 



Storzer & 
member. 



Hollenback), so it might be a very young cluster 



S-l M H 



|™-mh-Jr II iti 




1.L 1 



i i i i 

00 

_l L 



r|i O- Hi JiOIP ppl tnail 



nio-ina 

L A(V)=S.fl 



J 5J 

li-ihEiL-i i i, r»- i 



Figure 7(b): Observed if bandpass spectra of Trapezium sources. 
The stronger H2 emission lines in 047-436 are marked, together 
with a telluric feature at 2.0 /im. 

but an extinction of A(V) = 10 is required to produce a flat 
F„ spectrum, which might be more appropriate. Hence, the 
mass is highly uncertain and it is possible that it is a very 
low mass star rather than a brown dwarf. 

Object 019-108 also has high extinction and shows only 
weak water absorption at H but its K bandpass spectrum 
shows no sign of emission by H2 or by hot dust. This suggests 
that it is a background star rather than a cluster member, 
as discussed in Section 3.1. However the evidence is not con- 
clusive: circumstellar matter could be entirely swept away 
by the photoionising UV radiation field in 1 Myr (eg. 



4 ANALYSIS 

We analyse the spectra by comparison with both the 
AMES-Dusty-1999 model spectra calculated by Allard et 
al. (2000,2001) and with local field dwarfs. Since the spectra 
are dominated by water absorption bands, a water index of 
some sort is the natural tool to use to derive spectroscopic 
effective temperatures (T e //). We do not fit the full 1.4- 
2.5 /im spectra to those of field dwarfs or giants, after the 
fashion of Luhman et al.( 1998a) at optical wavelengths since 
the H bandpass profiles are clearly different (see Section 1). 
The synthetic spectra are believed to be most accurate at 
A < 1.7 /im (see Section 4.1). In addition, the most strongly 
temperature sensitive broad band feature at the observed 
wavelengths is the wing of the 1.35 /im water band in the 
1.4-1.7 fim region, the other bands being nearly saturated 
for all late M and L spectral types. Hence we define a water 
index W =F„(1.50 /an) /F„ (1.682 /im). The model spectra 
are rebinned to the same dispersion as the data and con- 
volved with the instrumental profile as measured in the Arc 
spectra. For the model spectra the 1.5 and 1.682 /im fluxes 
are averaged over 0.02 /im (or 9 data points) at both wave- 
lengths, while the observed spectra have their 'continuum' 
fitted with high order cubic splines to reduce the measure- 
ment errors. The 1.50-1.682 /im range covers most of the 
depth of this water absorption band, the data being less 
precise at shorter wavelengths. The model spectra peak at 
1.682 /im in the 2500-3300 K range and the dereddened data 
also peak at or very close to this wavelength. 

The choice of water index is always somewhat arbi- 
trary, so we also compute a reddening independent index 
of the type used by Wilking, Greene & Meyer at K band, 
but applied here to the H band spectra. This index is de- 
fined by Q = (F1/F2)(F3/F2) 1219 , where Fl is the F„ flux 
at 1.57 /mi, F2 is the F„ flux at 1.682 /im and F3 is the 
F„ flux at 1.79 /im. The 1.57 and 1.79 /tm wavelengths are 
chosen because they must be roughly equidistant from the 
flux peak for good precision and because the data is often 
of poor quality at A> 1.8/tm, and therefore not used in the 
continuum fitting process. We also tested an index similar 
to that of Delfosse et al.(1999), based on the slope in the 
1.51-1.57 /tm region. However, the measurement errors for 
the slope in this narrow region are generally too large for 
that index to be very useful here and K bandpass data is 
required to normalise the index as defined by Delfosse et al. 
The measured W and Q indices are listed in Table 3, along 
with spectroscopic temperatures derived by fitting the W 
indices to those of the models. Temperatures are rounded 
to the nearest 50 K, since the model grid has an interval of 
100 K. 

We assign approximate spectral types to each source 
(see Table 3), by the comparison of their W and Q indices 
with those of the local field dwarfs observed by Leggett et 
al.(2001), whose spectra types are defined using the 2MASS 
Kirkpatrick et al.(1999) system. The reddening independent 
Q index should be the more reliable for measuring water 
depth and for spectral typing but the W index generally 
produces very similar spectral types, as shown in the Table. 
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Figure 8(a): Dereddened spectra for Trapezium sources observed 
at both H and K. Noise in the 1.8-2 /j,m region is due to time 
variable telluric water absorption. 
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Figure 8(b): Dereddened spectra for Trapezium sources observed 
at both H and K. 

It is important to note that these are not real spectral 
types as might be defined at optical wavelengths - they are 
merely the spectral types that the H bandpass water depth 
would correspond to if the the Trapezium sources were high 
gravity field dwarfs. Low gravity sources have stronger wa- 
ter absorption for a given T e ff than field dwarfs. Hence 
the spectral types appear rather late for the corresponding 
spectroscopic temperatures in Table 3 and we would expect 
many of the L-type sources to have late M-type spectra at 
optical wavelengths. The behaviour of the W and Q indices 
as a function of spectral type is shown in Figure 9 for the 
field dwarfs. Most of the 15 field sources follow a linear rela- 
tion quite well but there is some interesting scatter (which 
is definitely real). Gdl65b is an L4 dwarf (in fact the pro- 
totypical object) but has very weak water absortion for its 



spectral type. The 2 M9 dwarfs in the Leggett al.(2001) 
sample show large scatter on either side of the best fit lines, 
particularly in the W index. Some of these differences might 
be attributed to dust effects: dust in cool metal rich atmo- 
spheres is expected to heat the photosphere; this may be in- 
fluenced by metallicity, rotation and age, as noted by Leggett 
et al. However, there is no clear difference between the water 
depths of Young Disk and Old Disk sources in the sample, 
categorised kinematically by Leggett et al.(2001). Another 
study by Leggett et al.(2000) showed little difference due to 
metallicity in the H bandpass when comparing Disk dwarfs 
and Halo subdwarfs but much weaker water absorption in 
Halo extreme subdwarfs, in which the H bandpass profile is 
dominated by H2 opacity in early M-type dwarfs. 

The spectral types of the Trapezium sources lie in the 
range M7-L5 for 16 of the 18 low mass sources (< 0.1 M@) 
which show strong water absorption. The 2 objects with 
later spectral types are faint and might have imperfect 
background subtraction (see Section 2). The fitted spectro- 
scopic temperatures (2100, 2350 K) are consistent with the 
new predictions for planetary mass sources by Baraffe et 
al. and by BM97 but are slightly low for these 2 sources 
with M~ 0.016 Mg at 1 Myr, unless they are younger, less 
massive objects of equivalent luminosity but greater size. It 
might be useful to define a new luminosity class or spectral 
type to describe these young, low gravity objects but since 
the gravities increase continuously over a period of several 
hundred Myr this is probably not a valuable concept. In- 
stead we regard these low mass sources as a low luminosity 
extension to the region occupied by T Tauri stars on the HR 
diagram. 

4.1 Synthetic Spectra 

The AMES-Dusty-1999 spectra appear to represent a signif- 
icant improvement over previous models (such as NextGen) 
of cool dwarfs at infrared wavelengths, as discussed by Al- 
lard, Hauschildt & Schwenke (2000, hereafter AHS). Models 
using the Miller- Tennyson list do not accurately reproduce 
the infrared profiles of either these data or field dwarfs (see 
Fig.9 of Leggett et al.2001). By contrast the AMES-Dusty- 
1999 spectra, which use the newer NASA Ames water line 
list, provide quite a good fit to the H bandpass profiles 
of the Orion data (see Figure 10(a-d)) and go some way 
toward reproducing to the plateau seen in spectra of field 
dwarfs. The synthetic spectra decline a little more steeply 
on the long side of the H bandpass than the data, unless 
T e // <;3300A". Such a high temperature would not be con- 
sistent either with the photometric temperatures obtained 
from the dereddened (1-3) colours in Paper I, or the evolu- 
tionary predictions of BM97 and Baraffe et al. As noted in 
Section 3.1, there appears to be incompleteness in the wa- 
ter line list for the absorption band centred near 1.95 /im. 
This compromises the predicted spectral profiles in the 1.7 
to 2.2 /im region and leads to small errors in the predicted 
H and K bandpass fluxes, as discussed by AHS. By con- 
trast the water absorption band centred at 1.35 /im is well 
fitted by the AMES-Dusty-1999 spectra (see AHS) and the 
fluxes and spectral profiles are believed to be more accu- 
rate at A < 1.7 /im. It is generally not possible to fit both 
sides of the H bandpass spectra simultaneously. This is one 
of the reasons why the short wavelength side is used to de- 
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Figure 9: W and Q spectral indices as a function of spectral type 
for Field Dwarfs. Best fit lines are plotted, excluding the divergent 
L4 source and the 2 M9-types. Including these would lead to later 
types for the coolest sources. Measurement errors are very small 
for these nearby sources. 

rive spectroscopic temperatures, the other being the greater 
temperature dependence of that part of the synthetic spec- 
tra. 

The discrepancy at 1.7-1.9 fira seen in Figure 10(a-d) 
could in principle be caused by hot circumstellar dust emit- 
ting detectable flux on the long wavelength side of the H 
bandpass. However, as shown in Figure 11, this would lead to 
strong and obvious excess emission in the K bandpass spec- 
tra, and some veiling of the water absorption near 2.5 /tm, 
which is not seen in most of the sources in Figure 8(a-b). 

The evolutionary calculations of BM97, D'Antona & 
Mazzitelli (1998) and Baraffe et al. predict surface gravi- 



ties between log(g)=3.5 and log(g)=4.0 for sources of brown 
dwarf and planetary mass in the 0.3-2 Myr age range 
expected for most Trapezium sources. AMES-Dusty-1999 
spectra in this gravity range reproduce the constant slope of 
the data at A < 1.68 /im and the wavelength of the flux peak. 
Higher gravity synthetic spectra show a flatter profile near 
the flux peak, and a steeper slope at shorter wavelengths. 
This is most easily seen in an Fa plot, shown in Figure 12, 
which may be compared with the composite Orion spectrum 
and field dwarf spectra in Figure 5(b). This is qualitatively 
more similar to the profile of field dwarfs (with log(g)« 5.5), 
although the abrupt change in slope observed in the 1.55- 
1.60 /im interval is not reproduced. The higher gravity mod- 
els also decline more steeply on the long wavelength side of 
the H bandpass than the low gravity ones, which leads to a 
greater discrepancy with the data at those wavelengths than 
is seen in the low gravity fits in Figure 10(a-d). 

In the K bandpass the behaviour of synthetic spectra is 
somewhat similar to that in the H bandpass but less clear 
cut. The low gravity models have a fairly constant slope on 
the wing of the water absorption band in the 2.05-2.25 ^m 
region, whereas high gravity models show a curved profile 
similar those in the H bandpass. The water vapour absorp- 
tion is less deep at K bandpass and the profile of the flux 
peak is apparently affected by H2 absorption longward of 
2.2 fj,m in some field L dwarfs (Tokunaga & Kobayashi 1999). 
One might argue that some of the spectra shown in Figure 
8(a-b), such as the high quality spectrum of 053-503, lack 
the flat profile seen in K bandpass F„ spectra of late M and 
L dwarfs (eg. Jones et al.1994) but there is more variation 
between sources than at H band, both in our data and in 
field dwarf spectra. Owing to these factors and the relatively 
poor quality of most of our K band spectra the water ab- 
sorption bands are not clearly distinguished from those of 
field dwarfs. 

The spectroscopic temperatures fitted in Table 3 gen- 
erally lie in the range 2150-3300 K for the 18/21 sources 
in the sample which have strong water absorption. All but 
2 sources have T e // > 2500 K, so there is insufficient dust 
in the photospheres to have a significant effect on the syn- 
thetic spectra. The fitted values of T e f f differ by some 450 K 
between log(g)=3.5 models and log(g)=4.0 models at the 
high temperature end, with the higher gravity models giv- 
ing lower temperatures. At 2500 K and below the same tem- 
peratures are fitted for both gravities. These temperatures 
are in good agreement with the evolutionary calculations of 
BM97. in the 0.008-0.080 M mass range, for which they 
predict 2150 < T eff < 3230 at 1 Myr, 2310 < T e// < 3290 
at 0.3 Myr and 2030 < T eff < 3180 at 2 Myr. The Baraffe et 
al. and D'Antona & Mazzitelli models generally predict very 
similar temperatures to BM97 (within 50-150 K) which in 
our view is less than the uncertainty in the calibration of the 
absolute temperature scale with observed spectral type. The 
log(g)=4.0 models produce better agreement for the higher 
mass sources of the sample than the log(g)=3.5 models. 

The spectroscopic temperatures fitted at log(g)=4.0 
also agree well with the photometric temperatures derived 
from the dereddened (I- J) colours from Paper I, by com- 
parison with the colours of field dwarfs (predicted for main 
sequence stars by BCAH using the NextGen model atmo- 
spheres). This is illustrated in Figure 13, which shows an 
excellent correlation between the two independently derived 
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Figure 10(a-d): AMES-Dusty-1999 temperature fits (solid lines) to 4 sources (dotted lines), fitted via the W index. The 1.4-1.7 /im 
region and the peak are well fitted. The 1.7-1.9 fim region cannot generally be fitted simultaneously, except in the warmest sources (see 

text). The discrepancy increases for higher gravities. 



temperature scales but a systematic offset toward lower tem- 
peratures of about 250 K in the photometric temperatures. 
The scatter is largely due to photometric measurement er- 
rors magnified by dereddening rather than the spectroscopic 
measurements, and fits in well with the estimated uncer- 
tainty of 0.25 mag in the dereddened (I-J) colours given in 
Paper I. The NextGen atmospheres are known to fit main 
sequence infrared colours quite well, (better then the AMES- 
Dusty models at T e f f > 2500 K) despite the very poor fit to 
the spectral profiles at H and K seen by Leggett et al.(2001). 
Hence the 250 K offset must be due either to the lower grav- 
ity of the sample or errors in the AMES-Dusty models. The 



NextGen spectra predict that (I-J) colours differ by only of 
order 0.1 mag between log(g)=3.5 models and log(g)=5.5 
sources, which corresponds to an interval in temperature of 
~50 K. However, the lower gravity also changes the pre- 
dicted (J-H) colours by ~0.1 mag, which leads to a larger 
difference in the (I-J) colours via the dereddening process. 
Hence, it is quite possible that the offset of 250 K is due 
solely to gravity effects, and the spectroscopic temperatures 
at log(g)=4.0 are correct. They are more consistent with 
the evolutionary calculations for age 1 Myr than the pho- 
tometric temperatures and are derived from the well-fitted 
1.4-1.7 /im region. Broad band (J-H) fluxes produced by 
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Figure 11: Hot dust. Strong emission by dust at 1200 K, cor- 
responding to excess E(H-K)=0.5 is added to an AMES-Dusty 
model. This can fill in the discrepancy at 1.7-1.9 /im but leads to 
a very obvious excess in the K bandpass. 

the NextGen and AMES-Dusty model atmospheres at these 
gravities are not yet more reliable than the assumption of 
main sequence colours, owing to the problems with the wa- 
ter line list referred to above. (The AMES-Dusty models 
produce most accurate colours for high gravity field dwarfs 
below 2500 K, while the NextGen models work best at higher 
stellar temperatures, above 3000 K). Hence, at present we 
continue to use empirical main sequence colours in the dered- 
dening procedure. 

There is a trend for less luminous, less massive sources 
to be cooler (see Figure 14), as indicated by the trend for wa- 
ter absorption to increase with decreasing mass. The BCAH 
models predict that water absorption increases monotoni- 
cally with decreasing temperature in these substellar pre- 
main sequence objects. A mass-temperature relation corre- 
lation is predicted by all the published pre-main sequence 
isochrones. The derived luminosities, and hence masses, are 
nearly independent of temperature owing to the very weak 
temperature dependence of the ( J-H) colour so this is not a 
trivial result. The correlation is imperfect, presumably be- 
cause of the dispersion in ages. 



4-1.1 Carbon Monoxide and other molecules 

The sythetic spectra also show a strong gravity dependence 
in the strength of the CO absorption bands near 2.3 /im, 
which is largest at low temperatures. This is illustrated in 
Figure 15, with K bandpass spectra of high and low gravity 
model atmospheres at 2500 K and 2900 K. We define a a 
CO index, K C o = (-F22 - F23) /F22 where F23 is the in- 
tegrated Fa flux in the 2.295-2.315 /im region and F22 is 
the integrated F x flux in the 2.2675-2.2875 /im region. The 
indices for models and data are listed in Table 3. All of the 
Trapezium sources have Kco < 0.12, i.e. an average absorp- 
tion depth across the first CO band of < 12%. The observed 
CO indices of 0.09-0.12 observed in 053-503, 014-413 and 
055-231 are in good agreement with the model predictions 
for log(g)=3.5 to 4.0 and 2500 < T eff < 3200 K. The data 
are much less consistent with high gravity models, particu- 
larly for the cooler sources. For example, we observe Kco = 
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Figure 12: Effect of gravity on AMES-Dusty models at H (F A 
spectra), (top) at 2500 K the [g]=3.5 spectrum (solid line) and 
[g]=4.0 spectrum (dotted line) are very similar but the [g]=5.5 
model (dashed line) has a flatter peak, (bottom) at 2900 K the 
[g=3.5] and [g=4.0] spectra begin to diverge. The [g=5.5] spectum 
requires a lower temperature, 2800 K, to produce similar water 
absorption and crosses the lower gravity spectra. 



0.119 ± 0.055 for 055-231, consistent with K C o = 0.12 at 
T ef f = 2575 K (with log(g)=3.5) but not K C o = 0.21 (with 
log(g)=5.5). Similarly, we observe Kco = 0.059 ± 0.035 for 
096-1944, roughly consistent K C o = 0.106 at T eff = 2700 K 
(log(g)=3.5) but not iv" C o=0.167 (log(g)=5.5). The model 
predictions for log(g)=3.5 and log(g)=4.0 are almost identi- 
cal. 

In 6/11 sources CO is not unambiguously detected in 
these low resolution spectra, since the drop in flux seen in all 
sources at A > 2.3 /im may be due to H2O as well CO, and 
the bandheads are not sufficiently resolved. We place upper 
limits on the Kco index for these sources in Table 3. Again 
these are consistent with low gravity models, but not with 
high gravity models in the case of the cooler sources, since 
CO depths of order 20% would be detected. This strong 
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Figure 13: Spectroscopic Temperatures [g=4.0], correlated with 
independent Photometric Temperatures for sources with mea- 
sured (dereddened) (I-J) colours. The correlation is excellent but 
the spectroscopic temperatures are higher by 250 K, as indicated 
by the dashed line for perfect agreement. 



absorption is observed in local dwarfs, as evidenced by our 
spectra of 3 red standards in Figure 16. 

In field dwarfs a fairly strong feature is often observed 
near 1.62 /tm, with uncertain identification at present. The 
AMES-Dusty-1999 spectra resolve this into a complex fea- 
ture with possible contributions from several lines of OH, 
Cal and Fel, with an OH doublet at 1.624 /im predicted to 
dominate in low gravity atmospheres. Giant spectra show 
the (v=6-3) CO bandhead at 1.619 /im in low resolution 
spectra, which increases in strength in very low gravity su- 
pergiants (Origlia, Moorwood & Oliva 1993) , but this is not 
predicted to be important at log(g) — 3.5 to 4.0. An absorp- 
tion feature is marginally detected at the 1.624 fim wave- 
length in the composite Orion spectrum in Figure 5(a-b) 
but it is not unambiguously distinguished from noise and is 
much weaker than the log(g) = 3.5 models indicate for the 
OH dominated features at T e// > 2800 K. At log(g) = 4.0 
the 1.624 fj,m feature is no stronger than at log(g) = 5.5, 
for the same temperature, so at that gravity the absence of 
a clear detection would be consistent with the data. (The 
AMES-Dusty-1999 grid has a sampling interval of 0.5 dex 
in gravity). 

Considering that the log(g) — 4.0 models also produce 
fitted spectroscopic temperatures more consistent with the 
evolutionary tracks than log(g) — 3.5 models, this might 
be taken as evidence that the gravities are closer to the 
higher value, though a finer grid in gravity would be needed 
to demonstrate this. On the BM97 1 Myr isochrone the 
3.5 < log(g) < 4.0 interval covers an age range from 
< 0.1 Myr to several Myr for 0.008 Mq planetary mass ob- 
jects, so a very precise gravity measurement would be needed 
to determine the age with useful precision. A slightly greater 
than average age might be expected for this spectroscopic 
sample of sources with low extinction. 
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Figure 14: Spectroscopic Temperature [g=4.0] vs. Mass. Masses 
are from the BM97 1 Myr isochrone. A weak correlation is ob- 
served and the temperatures are consistent with the essentially 
independent BM97 prediction. 

4-1.2 Atomic Lines 

The absence of the Nal line at 2.21 /im in 053-503 (our 
best K band spectrum) is also indicative of a low grav- 
ity atmosphere, as shown by the AMES-Dusty-1999 spec- 
tra in Figure 15. 053-503 shows only a weak local minimum 
at the 1-er level: this cannot be distinguished from noise 
but is consistent with the model. The weakness of Nal at 
2.21 /xm has been observed in a low gravity YSO by Luhman 
et al.( 1998b). In addition, the weakness of optical Nal and 
KI features in very young brown dwarfs has been observed 
by Bejar et al.(1999) in the cr-Orionis cluster; by Luhman, 
Liebert & Rieke (1997) in a brown dwarf in the p Ophiuchus 
region; and by Briceno et al.(1998) in pre-main sequence 
stars and brown dwarfs in Taurus. 



4.2 Comparison with cr-Orionis Spectra and 
IC348 data 

We note that the H and K bandpass spectra of 2 very 
low mass candidates in the cr-Orionis cluster (estimated age 
~5 Myr), plotted as low signal to noise smoothed spectra 
by Zapatero-Osorio et al.(2000), do not resemble those pre- 
sented in Figure 4(a-c). At H they show the plateau seen in 
field dwarfs and they also appear to be redder in the ratio 
of K : H bandpass fluxes than all of the sources plotted in 
Figure 8(a-b) with strong water absorption, except 014-413, 
which appeared to have an infrared excess due to hot cir- 
cumstellar dust. After converting the spectra to F„ the ap- 
parent ratios of flux near the bandpass peaks are K : H ~ 
1.6 and 1.2 for S Ori 47 and S Ori 60 respectively. The 
presence of the plateau feature might be taken as evidence 
that these 2 cr-Orionis objects are field L dwarfs, given that 
Zapatero-Osorio et al. expect ~33% contamination by these 
in their much larger survey area. However, there is sepa- 
rate evidence for the relative youth of 1 of the 2, S Ori 47, 
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Figure 15: Gravity effects on CO and Nal (F„ spectra), (left) 2500 K; (right) 2900 K. Both CO and Nal are much weaker in the low 
gravity [g=3.5] AMES-Dusty spectra (offset upper plots) than the [g=5.5] atmospheres. [g=3.5] and [g=4.0] are nearly identical at 

these wavelengths. 
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Figure 16: Field Dwarf spectra and 053-503. Late M and L type spectra of high gravity field dwarfs show very strong CO. 053-503 has 
similar CO to the M6 field dwarf but shows no Nal and stronger water absorption. 
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via low gravity features in its optical spectrum (Zapatero- 
Osorio et al. 1999) similar to those seen in Pleiades dwarfs. 
It would have to be an unusually young (~100 Myr) field 
L dwarf to have those optical features. In our view it is 
more likely that the H bandpass spectra and larger K : H 
ratios reflect the greater age of the <r-Orionis cluster than 
the Trapezium cluster. At 5 Myr, sources at or below the 
deuterium burning threshold are predicted to be 300 K or 
500 K cooler than at 1 Myr (Baraffe et al;BM97 calcula- 
tions respectively), corresponding to an effective tempera- 
ture (T e //) of 1600-1800 K for an object with mass 8 M. Jup 
and ~ 2200 K at 15 Mj up according to both calculations. 
This is substantially cooler than the Trapezium sources and 
similar to the temperature of L type field dwarfs which ex- 
hibit the plateau feature. This is well into the T e / / < 2500 K 
region in which dust in the photosphere starts to have a 
large influence upon the emergent spectrum, according to 
the AMES-Dusty-1999 models. The larger ratio of K : H 
band flux is also predicted for a lower temperature, being 
caused by the underlying Planckian shift of flux to longer 
wavelengths rather than circumstellar dust. There is a some 
indication of a plateau feature in the AMES-Dusty models 
at 2200 K in the 1.64-1.68 (im region (not shown) but it is 
much less broad than that observed by (Zapatero-Osorio et 
al.(2000) and weakens at lower temperatures. This might be 
due either to difficulties with the treatment of water or dust 
in the model or differences in the stratification of dust at 
this very early evolutionary stage. We emphasise that our 
H and K bandpass spectra of 3 field dwarf standards in 
Figure 5 and Figure 16 are in excellent agreement with the 
prior observations of Leggett et al.(2001), Jones et al.(1994), 
Tokunaga & Kobayashi and Delfosse et al. so we have no 
reason to doubt the reliability of the Trapezium spectra. 
Scattering in the circumstellar environment might affect the 
spectral profiles in subtle ways but would not erase features 
at specific wavelengths such as the H bandpass plateau. 



Najita, Tiede & Carr (2000) have discovered 6 sources 
at or below the deuterium threshold in the very young clus- 
ter IC348 (typical age 3 Myr), extending the work of Luh- 
man et al. (1998a). Although they did not have spectra, they 
were able to obtain approximate spectral types via narrow 
band imaging across the 1.9 /xm water band with NICMOS 
on the Hubble Space Telescope. They came to very similar 
conclusions to those presented here: the water absorption 
strength implied late M or L spectral types and they also 
argue that their sources lie above the main sequence and are 
therefore not background stars. This foreshadows our argu- 
ment in the next section. An interesting difference is that 
according to the AMES-Dusty models their 1.9 fim water 
index is stronger in low gravity pre-main sequence objects 
than in high gravity field stars of the same temperature. This 
apparent contrast with our earlier statement is due to the 
different broad band profiles of low and high gravity objects. 
The behaviour of water absorption strength as a function of 
gravity depends on the wavelength at which the water ab- 
sorption is measured. We will have to wait for the NICMOS 
Cooling System to be installed before measurements such as 
those of Najita et al. can be taken elsewhere, since telluric 
water absorption generally precludes good quality ground 
based observation at these wavelengths. 



5 CLUSTER MEMBERSHIP 

The reality of the population below 0.02 Mq has been 
called into question by Hillenbrand & Carpenter (2000), who 
suggest that lightly reddened background field dwarfs are 
present in the data, and are not distinguished from the pro- 
posed planetary mass population. They estimated that the 
total contamination by non-cluster members is 5% to K=18, 
using a modified Wainscoat model (Wainscoat et al. 1992) 
of the entire galactic stellar population. This type of calcula- 
tion is difficult, owing to substantial uncetainties in: (a) the 
infrared extinction on lines of sight through the cluster; (b) 
the scale height of the galactic M and L dwarf populations. 
This is relevant because the Trapezium lies « 122 pc above 
the plane, relative to the sun, which is more than 1 scale 
height above the Young Disk population but probably not 
the Old Disk. 

To address this issue we first look at the spectroscopic 
data to see whether they are consistent with a field dwarf 
population. Five dereddened H bandpass spectra of sources 
with 1 Myr masses M< 0.016 M Q and < A(V) < 2.3 
are shown in Figure 4(c). Although the signal to noise ra- 
tio (SNR) in these spectra is low, it is clear that all pos- 
sess deep water absorption on either side of the peak near 
1.68 /im, which would correspond to L-type spectra if they 
were main sequence dwarfs (see Table 3) . The spectral types 
are uncertain by a few subtypes due to the low SNR and the 
possibility of imperfect background subtraction but they are 
derived by a reddening independent method, so extinction 
does not affect the typing. Alternative water-based deriva- 
tions of spectral type could lead to earlier classifications, e.g. 
the Delfosse slope-based water index (see Section 4) which 
samples the steepest part of the water absorption in field 
dwarfs and not the plateau. However, we are confident that 
the H bandpass spectra of the very low mass candidates 
are of type M7 or later by any logical definition. For ex- 
ample, spectra of M0-M6 dwarfs do not decline significantly 
in flux/pixel (i.e. F\) on the short wavelength side of the H 
bandpass, going from 1.68-1.5 fim (see eg. Leggett et al.2000, 
1996), whereas the very low mass candidates all do, as illus- 
trated in Figure 17 for the candidates with weakest water 
absorption. Early M-type spectra rise from 1.7 to 1.25 fim, 
since water absorption is weak at these types in Disk popu- 
lation dwarfs and even weaker in Halo extreme subdwarfs. 

On the CIT photometric system an M7V dwarf has 
an absolute magnitude Mh = 10.22, according to a fit 
to sources with known parallaxes by Leggett (1992); or 
Mh = 10.26, averaging 2 sources from Leggett et al.(2001) 
with measured parallaxes in the Yale Catalogue (van Al- 
tena et al. 1995). Adopting the brighter value, this leads 
to m H = 18.44 for a source just behind OMC-1 at 440 pc, 
or rem — 18.57 on the UFTI magnitude system. This is 
more than a magnitude less luminous than the 3 very low 
mass brown dwarfs in Figure 4(c) and 0.6 to 0.9 magni- 
tudes less luminous than the 2 planetary mass candidates 
(with H dr =17.71, 17.95 on the UFTI system). Uncertainties 
in the dereddened H band fluxes for these 5 sources are ex- 
pected to be «0.2 magnitudes via the random photometric 
error and 0.15-0.2 magnitudes possible systematic error due 
to the uncertainty in the (J-H) colour of the isochrone. Later 
type dwarfs decline rapidly in brightness, eg. Mh = 10.90 
at M9, and 10.98 at L0 on the CIT system, as fitted by 
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Figure 17: F A spectra of 061-401 (0.011 Mq) and 038-627 
(0.014 M Q ) compared with field dwarf BRI0021-0214. The field 
dwarf spectrum is given added noise and smoothed for compari- 
son. 

Kirkpatrick et al.(2000) using sources with measured paral- 
laxes. Hence, our very low mass sources are almost certainly 
too luminous to be background main sequence dwarfs, even 
assuming zero extinction by the OMC-1 cloud. In fact, the 
backdrop of the OMC-1 cloud provides substantial extinc- 
tion in the region of the spectroscopic sample, as discussed 
below, which further reinforces this result. The luminosities 
are higher because they are physically much larger and have 
higher effective temperatures than field dwarfs with equiva- 
lent water absorption depths. 

The profiles shown in Figure 4(c) also resemble those 
of higher mass Trapezium objects rather than that of a late 
M or L type field dwarf, illustrated in Figure 17 with noise 
added to aid the comparison. The field dwarf has a plateau 
from 1.57-1.70 /xm, whereas the very low mass candidates 
appear to peak near 1.68 fim like their brighter counterparts. 

If background contamination were significant we would 
expect a range of lightly reddened M0-M6 dwarfs to appear 



in our sample, since these are the commonest stars which 
are luminous enough to be detected at d> 400pc, assum- 
ing very low extinction by OMC-1 in the outer regions of 
our survey. The earlier M-types are less common than the 
later types but would be detectable to a greater distance, so 
the distribution would be fairly even across spectral types 
M0-M6. Since they are not seen among our planetary mass 
candidates we suspect that HC have underestimated the ex- 
tinction for 2 reasons. The first reason is that the extinctions 
derived from our (J-H) colours show that the A(V) column 
extends beyond the spatial pattern indicated by C ls O maps, 
which are only sensitive to high density gas. High extinc- 
tion (up to A(V) = 14) is measured toward many sources at 
the edges of our survey (see Figure 18), which we attribute 
mainly to lower density molecular gas and also to the neu- 
tral cluster medium (O'Dell & Yusef-Zadeh 2000). The near 
infrared extinction measurements are direct, whereas values 
derived from radio measurements are always uncertain by a 
factor of a few. Nevertheless, there is radio waveband con- 
firmation of the extinction. The 13 CO and CS(l-O) maps 
of Tatematsu et al.(1993) sample the less dense molecular 
gas and show that the OMC-1 cloud has a low density wing 
extending far to the west of the main ridge in the cluster 
centre, beyond the location of our sources. The column den- 
sity of the wing corresponds to A(V)« 4 in the outer parts 
of the wing (using the transformations of Lada et al. 1994) 
for 13 CO and may be higher in the inner part where our 
sources lie: this region is not well resolved spatially from the 
very high density ridge in the map of Tatematsu et al. In ad- 
dition, the cold dust is seen directly in the 850 /im SCUBA 
data of Johnstone & Bally (1999), who measure continuum 
emission of 0.5-2.0 Jy/beam in the spectroscopic survey re- 
gion. Applying their adopted emissivity index (/3=2), tem- 
perature (20 K) and 100:1 ratio of gas:dust, we calculate 2- 
8xl0 29 Kg per beam (14 arcsec diameter). Adopting the re- 
lation A(V)/N(H)=5.3xlO -26 m 2 leads to 4.8 < A(V) < 19 
in the survey region, where the very low mass sources are 
detected. The uncertainties in grain emissivity, temperature 
and A(V)/N(H) mean that this is only an order of magni- 
tude estimate (eg. Ward- Thompson et al.1994) but it is con- 
sistent with the near infrared and 13 CO data. The second 
reason to expect significant extinction is that HC used an 
opticakinfrared extinction conversion based on the interstel- 
lar medium, which significantly underestimates the infrared 
extinction in Orion (Cardelli et al.) because of the presence 
of slightly larger dust grains in the environment of the molec- 
ular cloud. The substantial contribution of these slightly 
larger grains to extinction can be seen from calculations us- 
ing Mie theory (eg. Figure 15 of Lucas & Roche 1998). In 
some very young sources the measured infrared extinction 
may be local to the young star. However, comparison of Fig- 
ure 18 and Figure 3 does not show any less extinction in the 
outer region of our survey than in the full survey sample, 
which indicates that there is strong extinction by the low 
density molecular gas throughout the survey region. More- 
over, the L bandpass data of Lada et al.(2000) show that 
very young, embedded protostars are strongly concentrated 
towards the high density ridge, so the case for significant ex- 
tinction by lower density gas in the outer cluster is clear. It 
therefore appears that any background stars of type M7 or 
later would appear more than a magnitude fainter than all 
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Figure 18: Colour magnitude plot for the region of the spectro- 
scopic sample. High reddening is observed toward many many 
sources, similar to Figure 3. The new Baraffe et al. isochrone is 
plotted at left. 

the low mass candidates in Figure 4(c) and even background 
binaries would not be bright enough to cause confusion. 

Our assumption in Paper I that the OMC-1 cloud pre- 
vented any background contamination in the survey region 
at J band may prove to be an overstatement however. In 
Section 3.2.1, we discussed the 2 candidates near the stel- 
lar/substellar boundary (019-108, 066-433) in our spectro- 
scopic sample which suffer high extinction (A(V)=8.6, 12.6 
respectively) and have weak water absorption, consistent 
with late K-type or very early M-type spectra. If they are 
embedded YSOs then we can say that there is no background 
contamination in our sample, which is highly encouraging for 
future work in this cluster. It is very possible that they are 
background field stars, since their fluxes are consistent with 
that interpretation and the extinction is broadly consistent 
with what we might expect from the radio data. This would 
also be encouraging, since the high extinction towards these 
2 sources would indicate that background contamination can 
be avoided by studying sources which have low extinction. 

In Paper I we had assumed that lightly reddened back- 
ground sources are not present because of the large near in- 
frared extinction of cluster members and because the cluster 
profile at infrared wavelengths is apparently identical to that 
at visible wavelengths. A large population of field sources 
would be expected to produce a broader profile at infrared 
wavelengths, since these sources are uniformly distributed 
and more easily seen at the edges of the cluster, where ex- 
tinction is lower. This argument was first used by Hillen- 
brand & Hartmann (1998) but not at such faint sensitivity 
limits. We apply it here using our I and H bandpass data, in 
Figure 19, plotting the number of detections in each filter as 
a function of distance from 9i Ori C. Inspection of the pro- 
file gives no hint of a background population since the ratio 
of H:I bandpass detections is fairly uniform and may even 
decrease slightly between radii of 200 and 300 arcsec. This 
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Figure 19: Trapezium Cluster profile n(r)dr as a function of dis- 
tance from 8i Ori C. Solid lines show / band detections, dashed 
lines show H band detections. The profile is incomplete beyond 
160 arcsec because the survey region is not circular, but extends 
to the west. 

demonstrates that any population of faint sources (H ■> 17 
to be undetected at I) must be small in order to have no 
observable impact on the cluster profile. This method is not 
precise enough to exclude the presence of a small background 
population and there are selection effects in Figure 19 be- 
cause the detectabilty of sources increases more rapidly at 
large radii in the / bandpass than at H, owing to the greater 
effect of the nebular background at I. Nevertheless, the ab- 
sence of any obvious background population in the profile 
at H suggests that contamination is likely to be small at J. 

Could the planetary mass objects be foreground dwarfs? 
About half of the 15 candidates have sufficient reddening, 
(J-H) >0.85 (UFTI system) or (J-H) >0.94 on the CIT 
system, for us to be confident that they are significantly 
reddened. Only field dwarfs of type L2 and later can reach 
these colours according to Kirkpatrick et al.(2000) (or L6 
and later according to the smaller sample of Leggett et al. 
2001). These provide very small contamination, since mid 
and late L-types are only detected at small distances, sam- 
pling a much smaller volume than late M-types mj = 20.0). 
Kirkpatrick et al.(2000) The expected number of foreground 
objects is only 5 or 6 in the full survey region (~ 1% contam- 
ination) according to both the Tinney et al.(1993) luminos- 
ity function and the 2MASS data (Reid et al.1999), mainly 
M dwarfs in magnitude range 15 < J < 19. We expect to 
detect < 1 foreground L dwarf in the full survey region of 
40 arcmin 2 , using the space density estimate of Zapatero- 
Osorio et al.(2000). Although the volume density of field L 
dwarfs is uncertain, the number cannot rise by more than a 
factor of 3, since Zapatero-Osorio et al. calculated 33% con- 
tamination in their much larger survey area at very similar 
galactic coordinates (contamination cannot exceed 100%). 
Planetary mass candidates (< 0.012 Mq) are detected in no 
more than 60% or our survey area, where the nebular sur- 
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face brightness is sufficiently low. Hence, we derive an upper 
limit of 1 or 2 late M and L type sources which could appear 
to be planetary mass objects. 



6 AGE AND MASS OF SOURCES 

The interpretation of very faint cluster members as plane- 
tary mass objects might be incorrect if they are far older 
(> 20 Myr) than the bulk of the higher mass cluster mem- 
bers. The luminosity of a 0.020 Mq object does not de- 
cline to that of the 6 candidates with M < 0.008 M© un- 
til an age of 27 Myr, according to the BM97 calculations. 
The age distribution of the population has been studied in 
detail by Hillenbrand (1997) over the wider Orion Nebula 
Cluster, within a radius of 20 arcmin, using 2 evolutionary 
tracks in the HR diagram. This analysis concluded that most 
star formation in Orion has occurred within the past 2 Myr 
and that the youngest sources are highly concentrated to- 
wards the cluster centre, which is the region covered by our 
relatively small scale survey. The dataset did not include 
planetary mass sources but indicated that the lowest mass 
and highest mass sources are most strongly concentrated 
toward the cluster centre, with less concentration by inter- 
mediate (solar mass) sources. Specifically, Hillenbrand calcu- 
lated that < 2% of the population has an age >20 Myr, using 
the D'Antona & Mazzitelli (1994) tracks, which agree well 
with more recent calculations (BM97; D'Antona & Mazz- 
itelli 1998, Baraffe et al.). The Swenson et al.(1994) tracks, 
which indicated lower temperatures for brown dwarfs than 
more recent calculations, put the proportion at ~ 5%, which 
gives some indication of the uncertainty in this type of cal- 
culation. Conservatively, we assume that the old (> 20 Myr) 
population is only a factor of 3 less concentrated within our 
survey region (radius ~4 arcmin) than the young popula- 
tion (see Figure 22 of Hillenbrand 1997). This yields a con- 
tamination of 0.6%-1.7% by such old objects in our sample, 
or 4 to 11 objects spread over all masses in the sample of 
~ 650 sources detected at both J and H (including bright 
saturated sources). The more recent assessment of the age 
distribution in the core of the cluster by HC gave a simi- 
lar result. Expressed as N(log(M)), the IMF in Orion has 
been demonstrated to have a broad maximum somewhere 
between 0.1 and 0.6 M Q (Paper I, HC, Luhman et al.2000) 
and decline slowly at substellar masses. In the 0.03-0.08 Mq 
interval a relation N(log(M)) oc M _1 is consistent with our 
data in Paper I. If this trend continues below our complete- 
ness limit to 0.013 Mq, then we would expect there to be 
~ 0.16 to 0.44 old brown dwarfs with sufficiently low mass 
(~ 0.013 to 0.03 M 0) to masquerade as a planetary mass 
object. Even if the substellar IMF were flat we would not 
expect there to be more than 1 such source, whereas there 
are ~ 15 planetary mass candidates. In the event of a flat 
substellar IMF it would be difficult to argue for an abrupt 
cut off above the deuterium threshold. 

We note that the status of 1 brown dwarf candidate, 
014-413 as a youthful Trapezium cluster member has been 
confirmed by direct imaging of a wind bow shock interaction 
with the cluster medium in Ha in a Hubble Space Telescope 
survey (Bally, O'Dell & McCaughrean 2000) which unfortu- 
nately overlaps only a small part of the spectroscopic survey 
region. Outflow activity is correlated with accretion rate in 



YSOs, which decreases with age. Hence, this is probably 
a younger than average cluster member, which is also indi- 
cated by the presence of a large K bandpass excess in Figure 
8(a). 

K bandpass excesses are a good indicator of youth, since 
they are most common at ages < lMyr. Is there evidence 
for K bandpass excesses in our planetary mass candidates ? 
6/15 of our planetary mass candidates were detected in the 
H and K bandpass data of HC, which only partially overlaps 
our survey area. These can be dereddened with the addition 
of our J bandpass data. The H bandpass fluxes measured by 
HC agree well with ours in 3/6 cases, the measurements dif- 
fering by < 0.13 magnitudes, see Table 5. We note that the 
3/6 sources with poor agreement all lie near the edge of the 
survey of HC, and are believed to suffer from substantially 
shorter exposure times (25 s or 50 s total), whereas the other 
3 are believed to have the full 100 s exposure time. The poor 
agreement is therefore probably due to low SNRs in the HC 
data, given that all these faint sources are near the sensi- 
tivity limits of both surveys. Alternatively, large amplitude 
variability may occur in a small minority of sources, which 
has been demonstrated at stellar mass by Carpenter, Hillen- 
brand & Skrutskie (2001). If we accept that the 3 sources 
with good agreement have well measured H bandpass fluxes, 
then there appear to be large K band excesses in 2/3 cases. 
It is risky to compare photometry from different surveys 
in this way but these colour excesses, E(H-K)~ 0.6 magni- 
tudes, are large enough to be fairly convincing despite the 
uncertainty of order 0.1 mag in the (J-H) colour of the pho- 
tospheres. If they are genuine then this indicates that these 
are very young sources, probably with ages of <1 Myr, and 
therefore are very good planetary mass candidates. 

Oasa et al.(1999) were cautious in interpreting very 
faint sources in the Chamaeleon I cloud as planetary mass 
objects, despite the detection of large K bandpass infrared 
excesses requiring the presence of a substantial mass of hot 
dust in the accretion disk. However they found that even 
at 10 Myr, their presumed upper age limit for sources with 
substantial quantities of hot dust, a planetary mass was still 
indicated for their faintest source. 

At present, there is no easy way to actually measure 
the ages of individual sources because luminosity and tem- 
perature are degenerate with respect to mass and age, and 
because of the need to deredden sources. In principle the 
tracks of more massive sources might be separated on an HR 
diagram (older, more massive sources should be smaller and 
hotter than planetary mass objects of equal luminosity). At 
present, however the uncertainty in temperature measure- 
ment and prediction at these very early ages is too large to 
make such a distinction (see Table 3). A small bias toward 
greater ages might be expected for the spectroscopic sam- 
ple of sources with low extinction (see Section 4.3). Precise 
measurements of surface gravity via high resolution spec- 
troscopy of absorption line profiles may permit this at some 
future date but this will be difficult owing to the weakness 
of atomic lines in low gravity atmospheres. A promising 
alternative would be to use well modelled lines of simple 
molecules such as CO, and perhaps the Kco index, because 
some molecular features increase in strength with decreasing 
gravity. Another option might be to search for deuterium, 
the defining attribute of these 'free floating planets'. Deu- 
terium is predicted to be absent in sources more 10 7 years 
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old, for masses M> O.O3M0. This is illustrated by Bejar et 
al., using the calculations of D'Antona & Mazzitelli (1997), 
and examined in more detail by Chabrier, Baraffe, Allard, 
& Hauschildt (2000). 



7 GEMINI DATA 

We have conducted a pencil beam search for sources in the 
Trapezium even fainter than those deetected in our UKIRT 
survey, using the Gemini North Telescope with the Adaptive 
Optics system Hokupa'a and the infrared camera QUIRC 
The data will be described more fully in a future paper but 
preliminary indications are that no such ultra low luminos- 
ity sources were detected in 5 QUIRC fields containing ~ 25 
stars and brown dwarfs. This tends to support our previous 
suggestion thet the Luminosity Function (LF) in the Trapez- 
ium turns down somewhere near 8 Mj up , but this is far from 
conclusive because of the very small area surveyed and the 
existence of small scale spatial variation in the IMF, as seen 
in Figure 2). The decline of the LF should be investigated 
further to see whether it is a feature of the IMF or merely of 
the mass-luminosity relation. (An attempt at a much larger 
companion survey with the NAOMI Adaptive Optics system 
at the William Herschel Telescope was prevented by severe 
weather.) There are theoretical arguments for expecting few 
sources below 10 M Jup (Silk 1977; Adams & Fatuzzo 1996) 
but the physical processes in crowded clusters and multi- 
ple star systems may circumvent the obstacles described in 
those papers. 



cluster members, since they are too luminous to be back- 
ground field dwarfs and the modest extinction seen in most 
cases demonstrates that they are not foreground dwarfs. 
We would expect many of the L-type sources to show M- 
type spectra at optical wavelengths, since water absorption 
is deeper in low gravity dwarfs than high gravity dwarfs, ac- 
cording to the AMES-Dusty-1999 spectra. The fitted tem- 
peratures are higher than those usually associated with L- 
type sources. The trend of CO absorption to be significantly 
weaker at low gravities, for effective temperatures in the 
range 2100-2900 K, is a potentially powerful tool for iden- 
tifying very young brown dwarfs, provided there is an in- 
dependent means of measuring the effective temperature. 
This trend has been observed before in low gravity giant 
and supergiant stars (eg. Origlia, Moorwood & Oliva) at 
both 2.3 /im and 1.62 /mi. 

While the ages of the individual sources cannot be pre- 
cisely measured, it is very unlikely the 15 photometric plan- 
etary mass candidates could be much more massive sources 
with ages > 20 Myr, since statistically we would expect < 1 
such source in our dataset. There is also some evidence for 
large K bandpass excesses in 2 of the planetary mass can- 
didates, obtained by combining our JH data with the HK 
data of Hillenbrand & Carpenter 2000; however this must 
be confirmed by a single deep JHK survey. 



8 CONCLUSIONS 

Low resolution H and K bandpass spectra of brown dwarfs 
and planetary mass candidates show strong water absorp- 
tion bands, equivalent in depth to those seen in L-type and 
late M-type spectra, which is entirely consistent with the 
predictions of evolutionary tracks (Baraffe et al., BM97, 
DM98) at ages near 1 Myr. The new evolutionary tracks pro- 
vide self consistent flux predictions, permitting luminosities 
and masses to to be derived without the need for bolomet- 
ric corrections. The results are very similar to our previous 
work. However, the H bandpass profiles shown here are very 
different from those of field dwarfs with strong water absorp- 
tion, which appears to be a clear signature of youth and low 
mass. Strong evidence for the expected low surface gravities 
is provided by fitting the spectra to the AMES-Dusty-1999 
models. The low gravity models are quite successful at re- 
producing (a) the H bandpass profiles; (b) the weakness of 
CO absorption in the K bandpass spectra; and (c) the ab- 
sence of Na I absorption at 2.21 /mi, though with stringent 
limits in only one source. There is an excellent correlation 
between independently derived spectroscopic temperatures 
and photometric temperatures, with a relatively small sys- 
tematic offset of 250 K. 

Spectra of 2 planetary mass candidates and 3 very low 
mass brown dwarfs (M < 0.016 Mq at 1 Myr) show very 
strong water absorption consistent with L-type spectra, and 
have lower spectroscopic temperatures than the higher mass 
brown dwarfs, which is also predicted by the evolution- 
ary tracks. The late spectral types confirm that these are 
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Table 1 - Log of Spectroscopic Observations 
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25 nm resolution, K signal reduced by cloud 
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014-413 


14.44 


1.3 


0.056 


900 


960 


019-354 




013-306 


15.61 


6,1 


0.056 


3584 


- 


016-309 




096-1944 


14.96 


0.8 


0.042 


2240 


1792 


084-1939 




084-1939 


15.15 


0.6 


0.034 


2240 


1792 


096-1944 




084-104 


16.89 


0.0 


0.016 


2700 




095-100 




016-319 


17.26 


1.9 


0.016 


3584 


_ 


013-306 




038-627 


17.71 


2.2 


0.014 


3584 








061-401 


17.71 


0.0 


0.011 


5376 




066-433 




023-115 


18.41 


2.1 


0.008 


5376 


5376 


019-108 


K spectrum corrupted by nebulosity 


017-410 


12.23 


1.2 


0.30 


900 


960 


014-413 


serendipitous in slit 


073-205 


15.94 


11.8 


0.080 


960 




055-231 


serendipitous in slit 


178-232 


12.17 


14.5 ? 


0.4 ? 


240 






proplyd, higher mass 


G1406 


6.57 


0.0 






24 




M6V dwarf 


BRI0021-0214 


11.85 


0.0 




960 


180 




M9.5V dwarf, 25 nm resolution at H . 


DBD0205-1159 


13.59 


0.0 




720 


720 




L7V dwarf 



Notes: 

(a) H magnitudes are in the UKIRT - UFTI system. 

(b) The quoted masses are derived from the BM97 tracks, assuming an age of 1 Myr. 



Table 2 - H 2 emission in 047-436 



Measured A (/^m) Line Equivalent Width (A) 

1.945 (2-1) S(5) 5 ± 2 

1.956 (1-0) S(3) 18 ± 2 

2.032 (1-0) S(2) 4 ± 2 

2.122 (1-0) S(l) 15 ± 2 

2.202 (3-2) S(3), (4-3) S(5) 4 ± 2 

2.223 (1-0) S(0) 5 ± 2 

2.249 (2-1) S(l) 4 ± 2 

2.345 (4-3) S(3)? 5 ± 2 

2.407 (1-0) Q(l) 15 ± 3 

2.423 (1-0) Q(3) 9 ± 3 

2.453 (1-0) Q(5) 9 ± 3 

2.474 (1-0) Q(6) 7 ± 3 
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Table 3 - Spectral Types and Temperatures 



Source" 


Water Indices 6 


TefT(K) c 


Sp.Type d 


Mass (M ) 


Source Notes 




W 


Q 


[g]=3.5 


[g]=4.0 


W 


Q 






095-100 


U.UU ( XI. U± J- 


U.UJO XI.UJ.O 






i n 


LI 


0.090 




019-108 


n QQQ _i_ nn& 


n qq« 4. ni c: 


4850 


4900 


<M1 


<M1 


0.090 


high A(V), v. young source or background star 


055-231 




n 47Q 4. ni c; 


2600 


2550 


L3 


L4 


0.080 




069-210 


n 794 -U ni 7 


n c:Q7 4. n^1 


3400 


2950 


M7 


M9 


0.074 




068-020 




n kkx. _i_ 091 


3450 


3000 


M7 


LI 


0.073 




053-503 


n fi7i -u 04n 




OUJU 




i n 


M7 


0.073 




091-019 


n 77fi 4- 09^ 
U. 1 1 u xi .u^o 


n fil c; 4- 094 


3550 


3200 


M5 


M8 


0.071 




019-354 


n 717 + nn 

"J . I _L 1 1 .III.) 


n c:qft _i_ oi 4 


3400 


2900 


M8 


M9 


0.067 




066-433 


n 04 c; 4- 0^1 

U.lf-tO I.UJ1 


QQ7 4- Ofifi 


4450 


4500 


<M1 


<M1 


0.066 


high A(V), v. young source or background star 


043-014 


u.uji i • yj'-t'-t 


n 5 4- o^S 


2700 


2650 


LI 


LI 


0.061 




047-436 


0.785 ±.019 


0.921 ±.047 










0.06 ? 


H2 emission source 


014-413 


0.672 ±.014 


0.615 ±.043 


3100 


2750 


M9 


M8 


0.056 




013-306 


0.618 ±.031 


0.594 ±.018 


2650 


2600 


L2 


M9 


0.056 




096-1944 


0.632 ±.031 


0.562 ±.022 


2700 


2650 


LI 


L0 


0.042 




084-1939 


0.581 ±.050 


0.455 ±.035 


2550 


2550 


L3 


Lo 


0.034 




084-104 


0.471 ±.066 


0.396 ±.091 


2350 


2350 


L8 


L8 


0.016 




016-319 


0.359 ±.094 


0.394 ±.110 


2150 


2100 


>L8 


L8 


0.016 


poor sky subtraction ? 


038-627 


0.601 ±.072 


0.453 ±.044 


2600 


2600 


L3 


L5 


0.014 




061-401 


0.552 ±.054 


0.480 ±.085 


2500 


2500 


L5 


LI 


0.011 




023-115 




0.582 ±.150 








L0 


0.008 




017-410 


0.851 ±.047 


0.617 ±.030 


3850 


3900 


~M0 


L2 


0.30 


higher mass, unusual profile 


073-205 


0.693 ±.280 


0.469 ±.082 


3300 


2800 


M9 


L5 


0.080 


uncertain W index due to high A(V) 


178-232 




0.834 ±.017 










0.4 ? 


higher mass, proplyd, A(V) uncertain 



Notes: 

(a) The 3 sources below the line are not in the primary sample, see Section 3.1. 

(b) The W, and Q water indices are defined in the text. The quoted errors are based on apparent noise and do not include uncertainties 
in dereddening or the background subtraction. 

(c) Effective Temperatures are derived from the best fit of the AMES-Dusty-1999 spectra to the 1.45-1.7 fira data at two different 
gravities. 

(d) Spectral types are derived from the W and Q indices, calibrated by field dwarfs from Leggett et al.(2001) 



Table 4 - CO v=2-0 detections 



Source 




CO Detection 


053-503 


0.095 ±0.005 


yes 


014-413 


0.097 ±0.030 


yes 


096-1944 


0.059 ±0.035 


yes 


055-231 


0.119 ±0.055 


yes 


068-020 


0.063 ±0.030 


? 


084-1939 


0.032 ±0.047 


110 


091-019 


0.023 ±0.053 


110 


047-436 


0.034 ±0.026 


110 


019-108 


0.039 ±0.080 


110 


069-210 


0.153 ±0.082 


110 


019-354 


0.058 ±0.040 


110 


DBD0205-1159 


0.196 ±0.025 


yes 


BRI0021-0214 


0.154 ±0.007 


yes 


G1406 


0.102 ±0.007 


yes 



Note: Sources with no CO detection typically have Kco =0.02-0.03 d ue to the wing of the H2O absorption band H2O centred at 2.7 (im. 
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Table 5 - K bandpass excesses in Planetary Mass Candidates 



Source 


J 


rra 

LR 




K 


A(V) 6 


E(H-K) C 


AH d 


K band excess 


115-106 


19.32 


18.41 


18.30 


17.11 


2.49 


0.72 


0.11 


yes 


076-147 


18.67 


17.94 


17.93 


16.98 


1.11 


0.49 


0.01 


yes 


084-322 


18.74 


17.91 


18.04 


17.40 


1.90 


-0.02 


0.13 


no 


066-229 


19.97 


18.84 


18.15 


17.59 


4.31 


? 


0.69 


? 


137-531 


18.89 


17.90 


18.32 


17.28 


3.38 


? 


0.42 


? 


092-531 


18.72 


17.91 


18.46 


17.29 


1.78 


? 


0.55 


? 



Notes: 

(a) Hlr and Hue are the H bandpass fluxes measured in Paper I and HC respectively. 

(b) Visual extinctions are derived from (J-H) colours, with uncertainty of at least 1 mag. 

(c) An intrinsic colour of (H — K) = 0.38 is assumed for the photospheres. 

(d) AH is the discrepancy between the Paper I and HC H bandpass fluxes. 
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